Scaling behavior of transverse momenta distributions in 
hadronic and nuclear collisions* 

M. Rybczynski 1 ", Z. Wlodarczyk* 

Institute of Physics, Jan Kochanowski University, Swietokrzyska 15, 25-406 

Kielce, Poland 

AND 

Grzegorz Wilk § 

National Centre for Nuclear Research, Hoza 69, 00-681 Warsaw, Poland 

CM 
di 

^ C| It has been recently noticed that transverse momenta (pr) distributions 

observed in high energy production processes exhibit remarkably universal 
scaling behavior. This is seen when they are in some suitable variable, re- 
placing the usual px- On the other hand, it is also known that transverse 
momenta distributions in general follow a power-like Tsallis distribution, 
rather than an exponential Boltzmann-Gibbs, with a (generally energy de- 
pendent) nonextensivity parameter q. We now show that it is possible to 
choose a suitable variable such that px distributions of particles produced 
in proton-proton interactions in a wide energy range can be fitted by the 
same Tsallis distribution (with the same, energy independent, value of the 
g-parameter). Similar scaling behavior in nucleus- nucleus collisions is also 
observed. The possible dynamical origin of the q parameter used in these 
fits will be discussed. 
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It is well known [U [2j [3] that px distributions in general follow a two- 
parameter power-like Tsallis distribution [4j characterized by some energy 
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dependent nonextensivity parameter q, which for q — >■ 1 becomes the usual 
one-parameter exponential Boltzmann-Gibbs (BG) distribution^: 



h q (pt) = C q 
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(where C q is a normalization constant )q The transverse momenta distri- 
butions are then supposed to bring information on the thermodynamical 
properties of the production process, on the temperature T parameter in the 
case of exponential BG distributions, or, additionally, also on some intrinsic 
fluctuations presented in such systems and described by the nonextensivity 
parameter q in Tsallis distributions. 

It was recently shown that transverse momenta {pt) distributions ob- 
served in high energy production processes exhibit a universal scaling be- 
havior when presented in suitable variables, for example: 
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Fig. 1. (Color online) p + p, (p + p) data for transverse momentum distributions 
for different energies [HI HH1 E] plotted by using as scaling variables p' T (left panel) 
and u (right panel). 



• The variable p' T [J]- It is defined by demanding that px at energy W 

should be connected with p' T at energy W' via p' T = px (W /W)^+z; 
it can be reproduced by a two-parameter Tsallis fit, h q (p' t ) , cf. Fig. 
Q] (left panel) ; 

1 In phenomenological fits the Tsallis distribution is sometimes used interchangeably 
with the old two parameter, purely phenomenological, power-like parametrization, 
the so called "Hagedorn formula" [5]. Here we shall not discuss this possibility. 

2 Actually, as shown in [6], h q (pr) can fit all available data on pr distributions from 
RHIC to LHC, i.e., up to pr 200 GeV/c. The fact that data behave in such a way 
that they can be fitted by a simple two-parameter Tsallis formula for such a broad 
Pt range is really a phenomenon awaiting a proper understanding. 
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• The variable u = u{pt) = (p T f_6. PT • ^ is discussed in [8j, cf. Fig. Q] 
(right panel) showing h q (u); here 6 is an energy dependent parameter, 
6 = 6(a). 

In both cases q is energy independent, i.e., one obtains a kind of g-scaling 
phenomenon. 

It must be stressed at this point that essentially all data observed in 
high energy production processes can be represented by a universal scaling 
distribution, ip(z), [12, 13JE In Fig. H we demonstrate that these curves 
can also be nicely reproduced by two-parameter Tsallis fits h q (z). 




Fig. 2. (Color online) Tsallis fits to z-scaling plots presented by Tokarev. Left panel: 
spectra of charged hadrons produced in pp collisions at energies y/s = 19 — 2360 
GeV [13 ; right panel: spectra of different particles produced at y/s = 200 GeV 
(from talk by Tokarev - these proceedings). 

In this way we have reached our main point. Both h q (p' T ) and h q {u) 
(not to mention h q (z)) depart from the usual connection of the nonexten- 
sive Tsallis distribution with a thermodynamical approach (justified and 
advocated in [14]). Instead, as was already mentioned in [8], we opt for 
some more dynamical source of h q . As shown there, our scaling seems to 
indicate close connections with the so called preferential attachment (where 
T = T(pt)), seen in scale- free networks and dynamics leading to it [8]. In 
the case of variable u, this means that the parameter b can be positive or 
negative, depending on circumstances. This in turn can result in growth or 

3 Where, in short: ip{z) = —ixs/ [(dN / differ i n \ J 1 Ed 3 er /dp 3 , with J being the corre- 
sponding Jacobian of transformations from variables {p z ,pr} to {z,n}; z = zqQT 1 
where zo = ^/sr/[m - {dN c h/drj\o) c \ and S7 -1 is the minimal resolution at which 
a constituent subprocess can be singled out of the inclusive reaction, is the 

transverse kinetic energy of the subprocess consumed on production of mi and m,2, 
dN c h/dij\o is the multiplicity density of charged particles at n — 0, c is a parameter 
interpreted as a specific heat of created medium and m is is an arbitrary constant 
(usually fixed at the value of nucleon mass). 
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decrease of q. In our case q = 1.172, higher than obtained so far by using 
the usual h q (pr)- The question then arises of whether it is possible to get 
hq-+x(u) ~ exp (— u/uq) (by changing the sign of the parameter b in the 
definition of the variable u). In Fig. [3] we show a comparison of data for p+p 
and ^4+^4 for such scaling. They are different because of different uq and dif- 
ferent normalizations used. However, both curves coincide when multiplied 
by uo/A and when plotted for u/uq (they become independent of ^4), cf. 
Fig. [U They correspond to, respectively, b(^) = —0.085 + 0.115(y / s) -0 ' 2 
for p + p collisions and b(y/s) = -0.052 - 0.0002(^/s) a7 for A + A collisions. 
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Fig. 3. (Color online) p + p, {p + p) (left panel) and central A + A (right panel) 
data (from, respectively, [HI EH HI] and [TS] \T(^ H3 US] ) for transverse momentum 
distributions for different energies plotted by using the scaling variable u. 
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Fig. 4. (Color online) The same as in Fig. [3] but multiplied by uq/A and plotted 
for u/uq. Notice that curves for pp and ^4^4 coincide (left panel) and can be fitted 
by the same formula (right panel). The observed departure from an exponent for 
large u is connected with the effect of limitation of phase space and " conditional 
probability" discussed in [2]. 

We close with some remarks concerning the variable u. As already said, 
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by switching to this variable we depart from the usual thermal-like de- 
scription of such processes because u contains some dynamical input, for 
example: 

• One can write u/uq = pr/T e ff where T e ff is an effective tempera- 
ture: T e ff = T + T V (p T ) with T = u ■ (p T ) and T v = -b ■ u ■ p T . 
Such T e jf could be related to the possible px transfer, additional to 
that resulting from a hard collision, perhaps proceeding by a kind of 
multiple scattering proces^j. 

• T e j f also occurs in a description of the growth of the so called complex 
free networks which can then be applied to hadronic production ([8], 
cf. also [20] and references therein). 

Lets elaborate some more on the second proposition. Here one asso- 
ciates pt with the number of links in the quark-gluonic network assumed 
to be formed in the hadronization process. In this case their actual original 
energy-momentum distributions would be of secondary importance, since, 
because of their mutual interactions, they connect to each other and this 
process of connection has its distinctive dynamical consequences. One can 
think of such a process in the following way: We start with some initial 
state consisting of a number no of already existing (qq) pairs (identified 
with vertices in the network) and we add to them, in each consecutive time 
step, another vertex (a new (qq) pair), which can have ko possible connec- 
tions (links in network language) to the old state. Quarks are dressed by 
interaction with surrounding gluons and therefore " excited" and each quark 
interacts with k other quarks (has k links). Assuming that the "excitation" 
of a quark is proportional to the number of links k (which is proportional 
to the number of gluons participating in "excitation", i.e., existing in the 
vicinity of a given quark), the chances to interact with a given quark grow 
with the number of links attached to it, k. The new links will be preferen- 
tially attached to quarks with k. This corresponds to building up a so called 
preferential network, which evolves due to the occurrence of new (qq) pairs 
from decaying gluons. Such networks always result in a power-like behavior 
of suitable variables, in our case in pr- 
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